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ABSTRACT 

We have performed ah initio neutrino radiation hydrodynamics simulations in three and two spatial dimen¬ 
sions (3D and 2D) of core-collapse supernovae from the same 15 Mq progenitor through 440 ms after core 
bounce. Both 3D and 2D models achieve explosions, however, the onset of explosion (shock revival) is delayed 
by ^100 ms in 3D relative to the 2D counterpart and the growth of the diagnostic explosion energy is slower. 
This is consistent with previously reported 3D simulations utilizing iron-core progenitors with dense mantles. 

In the ^100 ms before the onset of explosion, diagnostics of neutrino heating and turbulent kinetic energy favor 
earlier explosion in 2D. During the delay, the angular scale of convective plumes reaching the shock surface 
grows and explosion in 3D is ultimately lead by a single, large-angle plume, giving the expanding shock a 
directional orientation not dissimilar from those imposed by axial symmetry in 2D simulations. We posit that 
shock revival and explosion in the 3D simulation may be delayed until sufficiently large plumes form, whereas 
such plumes form more rapidly in 2D, permitting earlier explosions. 

Subject headings: neutrinos — stars: evolution — stars: massive — supemovae: general 


1. INTRODUCTION 

That massive stars explo de at the end of their lives is well 
established observationally ( |Smartt|2009| ). Numerical simula¬ 
tion of core-collapse supernovae (CC SNe) has been less con¬ 
sistently successful than Nature (e.g., |Janka||2012[ |Hix et al. 
|2Q14t |Mezzacappa et al. ||20 1 5| ) . The ultimate source of the 
neutrino-driven explosion mechanism is the conversion of the 
gravitational binding energy of the core, collapsed to a proto¬ 
neutron star (proto-NS), and of matter accreted onto the proto- 
NS, into neutrinos that heat material behind the shock — re¬ 
viving it and expelling the stellar envelope as a supernova. 
This process is decidedly non-spherical. Neutrino heating 
above the proto-NS drives neutrino-dri ven convection (z/DC; 
|Bethe||1990t [Herant et al.||1992| |1994| ). Also excited are the 
low-or der modes of the sta nding accretion shock instability 
(SASI; Blondin et al.||2003| ). Asphericities in the shock sur¬ 
face channel the continuing accretion into distinct streams. 
All of these emergent behaviors are manifestly different with 
imposed axisymmetry (2D) than without (3D); therefore, we 
should expect 3D modeling to impact the initiation and sub¬ 
sequent development of CCSNe. 

Eully capturing the complex behaviors of the CCSN cen¬ 
tral engine numerically requires spectral neutrino transport, as 
the neutrinos are not in equilibrium with the fiuid and heating 
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is neutrino energy dependent, coupled to self-gravitating hy¬ 
drodynamics; i.e., spectral neutrino radiation hydrodynamics 
(zzRHD). Lew self-consistent, spatiall y 3D, spectral-uRHD 
CCSN simulations have been reported. Takiwaki et al. | ( |2Q 
computed a low-resolution 3D simulation of an 11 2-Mq pro 
genitor through the start of explosion. A 2. 7° angular reso 


lution follow-up using the same progenitor ( [Takiwaki et al. 
2Q14| ) also initiated explosions that subsequentl y developed 


more s lowly than corresponding 2D simulations. |Hanke et ah] 
(|2013|) computed a non-exploding 2.0° 3D simulation of a 27- 


Mq progenitor, though the same progenitor and physics did 
explode in 2D. 3D simulations of 11.2- and 20-Mq progen¬ 
itors with the same code and configur ation also failed to ex- 
plode ( [Tarnborra et al.|[2QT^ , though [Melson et ah] ( |2015a| ) 


found a delayed 3D explosion for a 20-Mq model with mod- 
ified opacities. This pat tern of delayed and failed explosions 
in 3D was suggest ed by|Hanke et al.|(|2012|) and later demon- 


strated explicitly ( |Couch||2013[ [Couch & Q’Connor ||20T4] ) 
using simulations with parameterized treatments of neutrino 


small differences or 

the opposite pattern dNordhaus et al. 

2010 

polence et al.||2013| 

|Handy et al. 2014 |Eemandez 


). Earlier work t 

)y Ery 

er & Warren (|2002| 

2004|) using 

non-spectral zzRHD found explosions in 2D and 3D to be sim- 


ilar. 

In this Letter we report on the first 440 ms of post-bounce 
evolution for a CCSN initiated from a 15-M0 progenitor 
in 3D and 2D using our multi-dimensional, zzRHD code 
Chimer;^ (Bruenn et al., in prep.) with modem neutrino 
interactions and general relativistic corrections to Newtonian 
self-gravity. We find that the shock revival occurs earlier by 
^100 ms in the 2D simulation relative to 3D, and that the 
growth of the diagnostic explosion energy is similarly accel¬ 
erated, potentially resulting in stronger explosions in 2D than 
3D. This is the first reported 3D explosion at 15 Mq, a rep¬ 
resentative mass often used for comparative studies, and the 
first for Chimera. In Sectionwe summarize our method- 


^ ChimeraSN.org 
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Figure 1. Mean (solid) shock radius for models C15-3D (green), C15-2D 
(black), and Cl5-ID (red) plotted versus time. Minima and maxima plotted 
with dashed lines. 

ology and initial conditions. An overview of the simulations 
is presented in Section with a focus on the differences be¬ 
tween the 2D and 3D simulations in Section [U We discuss 
our results in context in Section [^followed by a summary in 
Section |6l 


2. NUMERICAL METHODS AND INPUTS 


Initial cond itions are taken from the 15 Mq pre-supernova 
progenitor of|Woosley & Heger| (|2007|). The inner region 
(10,700 km; 2.32 Mq) is remapped onto 540 radial shells on 
logarithmic radial grid {dr /r) modified to track density gra¬ 
dients. Multi-dimensional simulations were initialized from 
a ID simulation at 1.3 ms after bounce by applying a 0.1% 
random density perturbation over radii 10-30 km, mimicking 
perturbations seen in simulations evolved through bounce in 
2D. The angular grid of the 3D simulation (C15-3D) was ini¬ 
tialized with a 180-zone (A(/) = 2°) 0-grid and a 180-zone 0- 
grid equally spaced in /i = cos i.e., equal solid angle. This 
0-gn& widens the pole-adjacent zones (A^ = i?sphA0sin6>) 
and therefore the time step . We evolve in spherical symmetry 
inside i?sph = 6 km until 45 ms after bounce (when prompt 
convection fades) thereafter setting i?sph = 8 km- With this 
grid, the pole-most zone is ?^8.5° wide resulting in a mini¬ 
mum leng th and time step ?^4 x larger than for a uniform 2° ^- 
grid (e.g., |Hanke et al.|2013 ). 300 ms after bounce, the O-grid 
was remapped in the 10 6^-zones closest to each pole (?^27°) 
to uniform spacing (A6> = 2.7°) and the ^-sweep at the pole 
was replaced by averaging, yielding similar time steps. The 
axisymmetric simulation (Cl5-2D) uses 270 uniform 6>-zones 
(A(9 = 2/3°). 


These are the third series of Chimera simulations (Series- 
C) and are substantially similar to the Series-B simulations 
dBruenn et al.]|MT^ |2^ hereafter |B2013| and |B2014| i. A 
more extensive desc ription of Chimera can be found in 
Bruenn et al.| ( |2014| ). The included microphysics are the 
same as for the Series-B models including the spherical GR 
terms in the gravity and transport. We solve the multi-group 
fiux-limited diffusion equations for all three fiavors of neu¬ 
trinos and anti-neutrinos with four coupled species: 

using 20 logarithmically 
spaced energy groups ae = 4-250 MeV, where a is the lapse 
function and e the comoving-frame group-center energy, in 
the ray-by-ray approximatio n. The neutrino-matte r interac- 
tions used are the full set of lB2014l We utilize the ILattiin^ 



Figure 2. Panel a: Diagnostic energy and Panel b: Mass of shocked 
cavity (solid), gain region (dashed ), and un bound region (dash-dotted) plotted 
in colors of Figure^ See text and |B2014| for definitions. 

& Swest^ ( |1991 ) EoS (incompressibility K = 220 MeV) for 


o >10^^ g cm ^ and an enhanced version of the [Cooper- 


stem 


g cm ^ 


_ , and in ou ter regions a 

14-species ct-network (|Hix & ThielemannJ 


(19851 EoS for p < 10 

119991 ) _ 

Relative to the Series-B simulations ( |B2013( |B2014| ), the 
neutrino transport solver now corrects for frame differences 
between shock-adjacent zones when computing the flux and 
fiux gradients (S. W. Bruenn et al., in prep.), permitting spher¬ 
ically symmetric Chimera simulati ons to track the late shock 
retreat of the reference simulation in [Lentz et al.| ( |2012] ). This 
improvement has a modest effect on the shock stalling radius. 

All times are given relative to core bounce. The proto-NS 
is defined as the volume where p > 10^^ g cm“^ and the 
shocked ‘cavity’ is the volume between the proto-NS and the 
shock. 


3. simulation OVERVIEW 

After remapping from ID, the multi-D simulations proceed 
in similar fashion: convectively unstable regions left behind 
by the shock progress through the Ee-core trigger prompt 
convection inside the proto-NS, similar to the axisymmetric 
Series-B simulations. 

Neutrino heating establishes a heating region extending in¬ 
ward from the shock to the gain surface, where net neutrino 
heating transitions to net cooling. Starting at ?^80 ms for both 
multi-D simulations, heating at the base of the gain region 
creates buoyantly unstable conditions, resulting in convective 
plumes rising against the continuing inflow. Rising plumes 
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Figure 3. a) Net neutrino heating in the gain region, b) (solid), Ve (dashed), and (dash-dotted) total luminosities at 1000 km. c) Neutrino heating 
efficiencies, d) (inward) Accretion rates at gain radius (solid) and shock (dash-dotted), e) Advection-heating time scale ratio, Tadv/^heat- f) Turbulent kinetic 
energy. Data for C15-2D is averaged with a 25-point boxcar (~8 ms). Plotted using colors of Figure[^ 


begin to affect the shock surface at ^95 ms for C15-3D and 
?^105 ms for C15-2D, as seen by the separation of the mini¬ 
mum and maximum shock radii (Figure[^ dashed lines). Over 
the next ^50 ms, both models become completely convec¬ 
tive within the shocked cavity. For C15-3D this results in 
a flat mean shock radius, i^shock, that rises gradually until 
^280 ms. For C15-2D, i^shock oscillates and grows faster, 
indicating earlier shock revival and explosion. The shock 
for Cl5-ID, which lacks multi-dimensional flows, reaches a 


maximum radius of «180 km at «80 ms and recedes there¬ 
after, typical of ID CCSN simulations. 

The shock in C15-2D expands rapidly from ^230 ms on¬ 
ward (Figure f^, with the diagnostic energ}p^F^+ (Figure]^) 
simultaneously becoming positive. E~^ surpasses 0.01 B oy 
250 ms and grows rapidly thereafter. For C15-3D, the first ev- 

following [ b 2014| is defined as the integral of the total energy (ther¬ 
mal, kinetic, and gravitational) in all zones of the cavity where locally posi¬ 
tive. 
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idence of potential explosion begins with an increased growth 
of i^shock at ^280 ms, accelerating after ?^350 ms, as the 
largest buoyant plume expands, leading to a small, but grow¬ 
ing . 

The explosion is clearly more energetic in C15-2D at all 
times (Figure [^). We evaluate the growth of over a com¬ 
mon period beginning when i^shock exceeds 500 km and end¬ 
ing 45 ms later. For C15-3D, i^shock passes 500 km at 393 ms 
when E^ is 0.034 B, which grows to 0.067 B at 438 ms when 
^shock is 735 km. For C15-2D, i^shock exceeds 500 km at 
278 ms when E^ is 0.041 B, which grows to 0.147 B at 
323 ms when i^shock reaches 900 km. Over this 45 ms com¬ 
parison period, the E^ growth rate is 0.73 B s“^ for C15-3D 
and 2.35 B s“^ for C15-2D, and the growth rate of i^shock 
in C15-2D is nearly double that of C15-3D. These indicate 
stronger growth of explosion in 2D. Further assessment of E^ 
growth requires significantly longer simulations: full satura¬ 
tion of E^ took ^1 second to reach ^1 B for the equivalent 
2D model (B15-WH07) of |B2014| 

The mass in the shocked cavity for C15-2D diverges from 
that of C15-3D at ?^150 ms (Figure |^; solid) and trends 
strongly upward by ^250 ms corresponding to rapid growth 
in i^shock and E^. Similarly, the mass in the gain region 
(dashed) grows from ^220 ms onward. For C15-3D, the turn¬ 
around in the mass curves is shallower, later, and occurs with 
less mass in the shocked cavity and gain region — all factors 
that correlate with weaker (or delayed) explosions. 


4. DIFFERENCES BETWEEN 2D AND 3D 


With successful explosions manifest in both 2D and 3D, 
the leading question is: “Why do axisymmetric models pro¬ 
ceed more rapidly (and more forcefully) to explosion than 3D 
counterparts?” 

From 50 ms to 150 ms, C15-3D shows larger total heating 
in the gain region (Figure |^), arising from greater heating 
efficiency (Figure |^; reheat; luminosity divided by net heat¬ 
ing rate) and neutrino luminosities ppigurej^), though mass 
accretion at the gain radius (Figure^; Mgain) is similar, re¬ 
sulting in larger i?shock than in Cl 5-^. At«150 ms, the ratio 
of advection and heating time scales (radv/^heau Figure [J), 
defined in |B2014[ grows past unity in bo th simulations, sig - 
naling the potential for thermal runaway ( [Buras et al.|2006| ). 
In C15-2D, Tadv/Theat gfows more rapidly, with large excur- 
sions driven by the oscillation of the shock along the pole. 
During this epoch, Mgain is larger in C15-2D, with large, 
positive excursions, correlated with favorable increases in the 
luminosities, r/heat. and Tadv/Theat. continuing through C15- 
2D shock revival (~250 ms). This favors earlier development 
of explosion in C15-2D, even though the luminosities remain 
higher in C15-3D. For C15-2D, luminosities and heating drop 
with accretion rate after shock revival at ^250 ms, while both 
remain noticeably higher in C15-3D. These measures of ac¬ 
cretion, luminosity, and heating are generally consistent with 
the early development of explosion in C15-2D. 

Recent work has examined the role turbulence can have in 


supporting the shock leading up to explosion (e.g., Murphy & 

" • - ' " " --- J2014' 


|Meakin|2011HCouch & Qtt|2015[|Abdikamalov et al.|20i'^ 

In Figure we show the kinetic energy associated with tur¬ 
bulence in the gain region. Turbulent kinetic energy 

1 


771A _ 

^turb — 


/ 

Jgain 


2P|v 


X|2 


dV. 


( 1 ) 


is defined for velocities, v^, integrated over the gain region. 
Lateral turbulent energy, computed by setting radial 


velocity Vr = 0 in defining Anisotropic turbulent en¬ 
ergy L^furb’ computed by removing the radial-shell mean 
(vr) from — (t’r)- The solid lines show the 

growth of L^lurb’ which begins growing at 100 ms, the onset 
of non-radial motions from convection, and continues after¬ 
ward. L’f^rb (dashed lines) is approximately four-fold larger 
in both simulations prior to shock revival. Both me asures are 
larger for C 15-2D than for C15-3D, consistent with |Couch &| 
Ott ( |2015| ), who posited that stronger turbulent pressure aids 


the development of explosions in 2D simulations. It is im 
portant to note that while L^^urb iiicasures the kinetic energy 
of disordered fiow, the relevant driver is convection driven by 
neutrino heating. This is especially important for where 
accretion streams and rising plumes are both deviations from 
the mean radial fiow at large scales. 

For multi-dimensional models (including those discussed 
above), there is a pre-explosion state with convective plumes 
rising to and distorting the shock. The fiow across the ac¬ 
cretion shock is diverted by shock geometry toward the local 
shock minima and then into accretion streams between the 
plumes. C15-3D, like previously reported 3D simulations, 
initially shows a large number of small plumes (Figure |^,b), 
whereas C15-2D (and most other 2D simulations) shows rapid 
d evelopm ent toward only a few large plumes (Figure]^). As 
in|B2014[ the primary polar plumes in C15-2D oscillate along 


the symmetry axis, in a manner consistent with the SASI, 
while neutrino-heated material continues to fiow into plumes 
from the bottom of the gain region, quickly triggering shock 
revival. In C15-3D, initially small rising plumes are pushed 
back by ram pressure at the shock and lack the persistence of 
the larger polar plumes in 2D models (see animated version 
of Figure!^). 

As C1>3D progresses toward explosion, the angular scale 
of the plumes slowly grows and, eventually, the largest of 
these plumes expands continuously and drives shock revival. 
The growth of plume scale can be seen at the surface (anima¬ 
tion of Figure]^) and deeper (150 km; Figure [^. The growth 
of the plume scales clearly precedes shock revival at ^350 ms, 
when the shock is clearly expanding and upfiow at 150 km is 
dominated by a single large pl ume with a few sm aller plumes. 
Parameterized simulations of |Femandez| (|20 1 5| ) show merg¬ 
ers of plumes and those of |Handy et al.| ( |2014| ) also found de¬ 
creasing plume numbers preceding, and after, shock revival. 
In 2D (including C15-2D), the axis-focused SASI and the ten¬ 
dency for convective vortices to merge ensures that the Targe- 
plume-state’ is reached earlier (by 150 ms in C15-2D). In 3D 
this state takes time to develop and appears to delay the ex¬ 
plosion. 

We have not yet completed a comprehensive analysis of the 
driver of plume-scale growth, but some observations from our 
simulation are revealing. As plumes grow in 3D simulations, 
the associated accretion streams are displaced farther from the 
central axis of the plume and are therefore less likely to in¬ 
terrupt the infiow of additional buoyant material, providing a 
‘self-shielding’ effect that permits further growth. Addition¬ 
ally, the curvature of the shock surface relative to the radial 
inflow reduces the effective ram pressure normal to the shock 
surface at locations farther from the plume’s axis and thus fa¬ 
vors lateral expansion of the plume. Together, these effects 
could drive the growth of larger plumes, which eventually are 
able to expand continuously, with material injected from the 
bottom of the gain region as in 2D simulations. 

The effectiveness of this mechanism may be seen in the 
growth of the plume in the upper-left comer of the animated 
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Figure 4. Specific entropy (ke baryon“^) at 200, 300, and 400 ms with 400-km scale bars in each panel. Column a (left): Volume rendering for C15-3D 
using a fixed transfer function, highlighting rising plumes. Column b (center): Polar slice through C15-3D, aligned with Column a. In upper two panels (200 and 
300 ms), the 180° (/)-shift between upper and lower halves is exaggerated by the 8.5° zone at the pole. 400-ms panel shows effect of transition to (^-averaging at 
pole. Column c (right): Entropy in a polar slice through C15-2D with color scale matching Column b at each epoch. Animated version of Column c available at 

ChimeraSN.org 
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Figure 5. 150-km slice of radial velocity in cm s ^, showing development of plumes. Upflows (downflows) are red (blue) saturated at +(-)10000 km s ^. 
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Figure]^, which precedes the onset of explosion in 3D. The 
growth of this plume (also visible in the 200 and 300 ms pan¬ 
els of Figure |%,b) eventually covers a significant portion of 
the shocked volume and shock surface solid angle. The lateral 
expansion of this plume diverts the accretion fiow of material 
striking it farther and farther from its axis, which contributes 
to the formation of a strong accretion region on the opposite 
side of the proto-NS at ^350 ms, clearly seen in the increased 
density accreting opposite the largest plume. The growth of 
the leading plume, and the stron g accretion opposite it, give 
C15-3D a preferred axis (see also Dolence et al.|2Q13) . 


5. DISCUSSION 


The role of dimensionality in CCSN modeling has been 
of intense interest recently starting with neutrino Tightbulb’ 
models that paramete rize heating with a fixed luminosity. 


Nordhaus et al.| ( |2Q10| ) found that 3D enhanced the potential 
for explosions, explaining the effect through increased dwell- 


timefor parcels in the neutrino heating region. Hanke et ah] 
( [2012| ) were unable to confirm that result and for their base 3"" 
simulations found little difference between 2D and 3D. They 
found that improved angular resolution enhanced explosion 
in 2D and inhibited explosions in 3D, which they attributed 
to the differences in the turbulent cascade and the action of 


the SASI. |Couch| ( 2Q13| ) found that 3D diminished the po¬ 
tential for explosions, while following the details of the first 
study carefully. Subsequent studies with other parameteriza- 
tions have mixed results. Parameterizations necessarily omit 
physics that might be critical to the nature of the CCSN mech¬ 
anism, thus the importance 2D versus 3D must be made rela¬ 
tive to simulations containing all needed physics. Most z/RHD 
simulations, including these, have shown that 3D diminishes 
the potential for explosion. 

The earlier revival of the shock for C15-2D is consistent 
with the other z^RHD simulations in the literature. Axisym¬ 
metric simulations of 11.2, 2 O-M 0 ( [Tamborra et al.||2013| 


2014 


Hanke 2Q14| ) and 27 -Mq progenitors ( [Hanke et al. 


2013 ) produce explos ions where their 3D counterparts do not. 


Melson et al.| ( |2015a| ) obtained a delayed explosion in 3D by 
altering neutrino-nucleon scattering for the 20-Mq progeni¬ 
tor. In the 3D 20- and 21-Mq models, they una mbiguously 


demonstrate the sp iral SASI mode identified by |Blondin & 


|Mezzacappa| ( |2007 , where as C15-3D resembles the neut rino 
dominated simula tions of [Abdikamalov et aL| ( |2014| ) and 
[Fernand^ ( |2015| ). In their multi-D simulations, the shock 
recede s, in the manner of I D simulations. Shock contraction 
in the |Hanke et al.| (|20131) 21 -Mq model ends, when, after 
accreting a composition interface that drops the (shock) ac¬ 
cretion rate from ^0.75 to «0.25 Mq the shock rapidly 
expands to near its previous peak. Following expansion, the 
shock revives in their 2D simulation, but contracts in their 
3D simulation. Likewise, in the 20-Mq simulation, a drop 
in accretion from «0.8 Mq to ?^0.4 Mq at 250 ms 
expandes i^shock without explosion, and in the II. 2 -M 0 sim¬ 
ulation a relatively fiat accretion of ^0 .2 Mq s~^, results in a 
gr adual decline of the 3D shock radius ( Tamborra et al.|2Qf4 ). 
In Melson et al. ( [2015a| ), the same accrkion drop in a 20-Mq 


3D simulation t riggers an explosion using enhanced neutrino 


heating. As in |B2014[ sudden decreases in accretion at the 
shock, Mgh, did not trigger shock revival for this progenitor in 
2D. For our models at shock revival, Mgh is ^0.7 Mq s~^ for 
C15-2D (^250 ms) and «0.5 Mq 8“^ for C15-3D (^350 ms) 

and smoothly declining (Figure [M; d ashed lines). _ 

In the multidimensional mod^s of |Takiwaki et al.| ( |2014| ). 


^shock rises smoothly beyond maximum ID i^shock, with the 
2D simulations growing faster than 3D, though these mod¬ 
els include less physics than the work of the Garching group 
discussed above and than our simulations. In their highest 
resolution simulations all of the 2D runs and the longest-run 
3D simulation show steepening of the i^shock curve, as seen 
in our models. _ 

In co ntrast to the above and to our results, [Melson et al.| 
( |2015b| ) find a modest increase in explosion energy for 3D 
versus 2D for a low-mass 9.6-Mq Fe-core progenitor for 
which they also obtain an explosion in ID. They attribute 
this difference to reduced cooling due to increased dissipa¬ 
tion of the accretion fiow above the gain surface. Given the 
progenitor differences (ID explosions, plumes not reaching 
the shock, etc.) these differences relative to our simulations 
are not surprising. 

For the remaining 3D z/RHD simulations, including C15- 
3D, it appears that there is some sensitivity to neutrino heat- 


of the Garching group ([Hanke et al.[ 

2013 

[Tamborra et al. 

2013 2014[ 

i [Melson et al.|2015a[) all e; 

xhibit declining T^shock 
y sharp declines in ac- 

that are onl 

y (temporarily) reversed b 


cretion. Only by modifying the neutrino opacities to account 
for strange quarks does the ^10-20% increase in neutrino 
heating turn the shock reversal into an explosion. Ostensibly, 
the Garching code is most similar to Chimera of all multi¬ 
dimensional z^RHD CCSN codes, though there are some dif¬ 
ferences: we do not use the ^ pai r-conversion 

process that increas es cooling ( jBuras et al.|2003b ; we use the 
|Cooperstem| ( |1985| ) nuclear EoS in the dynamically important 
shocked cavity; they use variable Eddington tensor transport 
along rays while we use flux-limited diffusion; and there are 
likely differences in numerical dissipation of fiows and turbu¬ 
lence from angular and radial resolution and radial grid mo¬ 
tion. The effects of these differences can only resolved by 
code-to-code comparison. 

6. SUMMARY AND CONCLUSIONS 

We have computed two identically initialized CCSN simu¬ 
lations, one in axisymmetry (Cl5-2D) and the other without 
imposed symmetry (C15-3D), for 440 ms after core bounce 
using the Chimera CCSN code. The shock revives fairly 
quickly (^250 ms after bounce) for C15-2D, following large 
£ = 1 oscillations of the shock from the SASI and buoy¬ 
ant plumes along the poles, similar t o previous 2D simula¬ 
tions with Chimera dB2013t|B2014|). Immediately preced¬ 
ing shock revival, C1 d- 2D exhibits rnore neutrino heating in 
the gain region and higher accretion at the gain surface than 
C15-3D. For C15-3D, the shock revival occurs ^100 ms later, 
while accretion and convection continue. Though these sim¬ 
ulations do not extend long enough to determine final ex¬ 
plosion energies, the energies reached in C15-3D, and their 
growth, make it likely that a lower final explosion energy will 
be reached in 3D, or it may take longer to reach a similar en¬ 
ergy relative to 2D. 

The shock is revived in C15-3D through the lead of a single, 
large-angle plume that results in expansion with a preferred 
axis, not wholly dissimilar to the axis-imposed structure of 
2D simulations. Based on examination of the buoyant plumes 
and their angular growth in C15-3D, we speculate that the 
development of fewer large-scale individual plumes may be 
necessary in 3D to provide the buoyancy needed to overcome 
the accretion ram pressure and relaunch the shock. 

We see evidence for larger turbulent kinetic energy con- 
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tributing to the earlier shock revival for C15-2D, but further 
examination is required. Though these are among the best- 
resolved spectral-z/RHD CCSN simulations reported in the 
literature, we can not be certain the accretion streams and tur¬ 
bulence are adequately resolved. We plan more extensive cov¬ 
erage of these issues, including an examination of resolution, 
in subsequent publications. 
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